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Abstract. We evaluate the total integral from negative infinity to positive infinity 
of all global solutions to the Painleve II equation on the real line. The method 
is based on the interplay between one of the equations of the associated Lax pair 
and the corresponding Riemann-Hilbert problem. In addition, we evaluate the total 
integral of a function related to a special solution to the Painleve V equation. As 
a corollary, we obtain short proofs of the computation of the constant terms of the 
limiting gap probabilities in the edge and the bulk of the Gaussian Orthogonal and 
Gaussian Symplectic Ensembles that were obtained recently in [3] and |T8]- We also 
evaluate the total integrals of certain polynomials of the Painleve functions and their 
derivatives. These polynomials are the densities of the first integrals of the modified 
Korteweg-de Vries equation. We discuss the relations of the formulae we have obtained 
to the classical trace formulae for the Dirac operator on the line. 



Mathematics Subject Classification: 33E17, 35Q15, 15A52 



Total integrals of global solutions to Painleve II 



2 



1. Introduction 

In this paper we compute the total integral, or integral from negative infinity to positive 
infinity, of all global solutions to the Painleve II equation on the real line (modulo an 
additve factor of 27riZ for one case; see Theorem 13.21) . If the solutions do not decay 
sufficiently fast as x — > ±oo then appropriate terms from the asymptotic expansion of 
the solution are subtracted off to make the integral convergent. One of the motivations is 
to give a new, short proof of the constant terms (first computed in [1]) in the asymptotic 
expansions of the distributions of the largest eigenvalue of a GOE or GSE matrix in the 
edge scaling limit. This employs the total integral of the special Hastings-McLeod 
solution (see Theorem 12. 2j) . In addition, in Section [5] we compute the total integral of a 
function related to a special solution of the Painleve V equation. This allows us to give 
a short proof of the constant terms (first computed in [18]) in the asymptotic expansions 
of the limiting gap probabilities in the bulk for a GOE or GSE matrix. In the last two 
sections we evaluate the total integrals of the polynomials of the Painleve functions and 
their derivatives that are produced by the densities of the first integrals of the modified 
Korteweg-de Vries equation. The evaluation of these integrals, although much simpler 
than the evaluation of the total integrals of the Painelve functions themselves, allows us 
to introduce the Painleve analogs of the classical trace formulae of the scattering theory 
(see equations (I252p in Section [7]). 

The homogenous Painleve II equation 



can be solved via a certain Riemann-Hilbert problem (see [22] ; also see [201 EH EH EH] 
for the derivation and for the history of the subject). Define the six rays '■ = 
{e 1 ( 2fe ~ 1 ) 7r ' 6 ]R + }, k = 1, . . . , 6 oriented outwards from in the complex plane. On each 
7^ define the jump matrix Sk as shown in figure [TJ The complex constants Si, s 2 , and 
s 3 satisfy 



u xx (x) = 2u 3 (x) + xu(x) 



(1) 



si - s 2 + s 3 + sis 2 s 3 = 0. 



(2) 




Figure 1. The Riemann-Hilbert problem for Painleve II. 
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Solving the Riemann-Hilbert problems means finding a 2x2 matrix valued function 
^(A; x) such that 

' \1/(A; x) is analytic for A ^ 7^, k — 1, ... ,6 
V I / + (A; x) and ^_(A; x) are continuous for A G 7^ k = 1, . . . ,6 
\I/+(A; x) = ^-(A; x)S k on 7 fc , with S k defined in figure m ( 3 ) 



*(A; x)e e(A;z)ff3 =I + o(jj as A 



00. 



Here ^+(A; x) and \I/_(A; x) denote the nontangential limits of \I/(A; x) from the left and 
right sides of the jump contour, respectively, and 

6(\;x) :=i(^\ 3 + x\^ . (4) 

If \1/ exists, 

u(x) := 2 lim (A*i 2 (A; x)e- e ^ x) ) (5) 

A^oo 

is a solution of (Q. Indeed, the Riemann-Hilbert problem is always uniquely solvable 
(the solution is a meromorphic function of x) and the map 

{si, s 2 , s 3 } — ► {set of all solutions of ([I])}, (6) 

defined by formula (jSJ), is a bijection (see Theorem 3.4, Theorem 4.2, and Corollary 4.4 
in [22]). 

We define the Pauli matrices as 

01 '= ( \ I ) ' a2 := ( i T ) ' " 3 := ( I -1 ) ■ (7) 

It is easy to check from the Riemann-Hilbert problem that ^(A; x) satisfies the Lax pair 
d 

= (-i(4A 2 + x + 2u 2 )a 3 - Au\a 2 - 2va 1 )^ (8) 

OA 

d 

—^ = (-iXa 3 -uo- 2 )^. (9) 

Here the function v(x) satisfies v(x) = u x {x) = dw(x)/dx = u'(x). The Painleve II 
equation (pQ) is indeed the compatability condition for this overdetermined syste 

Let ^/fc(A;x) indicate the function ty(\;x) restricted to A G Q k , where the regions 
Qk are defined in figure [TJ The x differential equation is particularly simple when 
A = 0: 

4-P(x) = -u(x)<r 2 P(x), P(x) := lim^ fc (A;x) (10) 
ax a^o 

f Note that here we use the Lax pair from p. 174 of [22 . This Lax pair differs from the original Lax 
pair suggested in [20] and reproduced on page 161 of [22j by a matrix conjugation with the matrix 
exp(if cr 3 ). 
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for some k, where A approaches in Q^. This limit is well defined since \l/fc(A;x) takes 
continuous boundary values. The general solution of ffTUl) is 



P(x) = e~ u ^P(a) = ( cosh f («^) ^U(a,x) ) p(a) 

-ismh U(a, x) cosh(7(a, x) 



[ID 



where a is a constant and 

U(a,x) := / u(y)dy. (12) 

J a 

Hence 



coshU(a,x) isinh?7(a, x) 
— i sinh U(a, x) cosh U (a, x) 



P(x)P(a)-\ (13) 



If we take x — > +oo and a — > — oo, this yields a relation between the total integral 
J^u^dy and the solution of the Riemann-Hilbert problem ([3]). Therefore, by 
analyzing the Riemann-Hilbert problem asymptotically as x — > ±oo using the Deift- 
Zhou steepest-descent method, we can compute the total integral. 

The asymptotic analysis asi-> ±oo for the Painleve II Riemann-Hilbert problem 
has been worked out in [H] and [22jij]. Most of the asymptotic analysis we will need in 
this paper is carried out in these references with the exception of the 0(x~ l ) term in the 
generic purely imaginary global solutions, which we compute in Section HI Nevertheless, 
the value of the solution at z = 0, P(x) := liniA^o ^k, has not been specifically addressed 
before, and in the subsequent sections we compute this term explicitly. We adopt 
the notation in [22] except when computing the total integral of the Hastings-McLeod 
solutions, when it is convenient to follow [II] where the original Riemann-Hilbert setting 
of [20] is used. We remark that the monodromy data (p,q,r), jump matrix Vnz, and 
solution rrv- l > to the Riemann-Hilbert problem in [H] are related to those in [22] by 

p=is 3 , q = is x , r = -is 2 , 

V DZ = e-iWV^Ve^e 1 * */ 4 , (14) 

m (l) = e -i7ro-3/4^ e 0o-3 e +i7ro- 3 /4^ 

Note that the phase factor e ecr3 appears in the normalization condition in [22] and in 
the jump matrices in [H] . 

We conclude the introduction with the following useful observation that relates the 
solution corresponding to monodromy data (s\, s%, s 3 ) to the solution corresponding to 
monodromy data (— s\, —S2, —S3) (cf. [33] and Chapter 11 of [22J). 

Lemma 1.1. If u(x; s%, S2, S3) is the solution to (QJ) with monodromy data (s%, S2, S3), 
then 

u(x; si, s 2 , s 3 ) = -u(x; -s u -s 2 , -s 3 ). (15) 

I We refer to the introduction of [55] for a detailed historic review on the asymptotic analysis of the 
Painleve equations via the Ricmann-Hilbcrt-isomonodromy method. 
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Proof. Define 

*(A; x) : = e imT3 / 2 $(A; x)e- [na3/2 . (16) 

Then \l/(A;x) satisfies the Riemann-Hilbert problem ([3]) with the jump condition 
replaced by \I/+(A; x) = \I/_(A; x)e m<T3 ' 2 Ske~ r * <T3 '' 2 on 7*.. The only effect this conjugation 
has is to change (si,s 2 ,s 3 ) to (— s 1; — s 2 , — s 3 ) in the jump matrices. Therefore, if 
\I/(A; x, Si, s 2 , S3) is the solution to the Riemann-Hilbert problem (J3j) with monodromy 
data (si,s 2 ,s 3 ), then \I/(A; 2, Si, s 2 , S3) = ^(A; x, — Si, — s 2 , — s 3 ) by the existence and 
uniqueness of the solution of the Riemann-Hilbert problem. From ([HJ), 

u(x; si, s 2 , s 3 ) = 2 lim (A$i 2 (A; x, -s u -s 2 , -s 3 )e~ e(x '' x) ) 

A^oo 

= -2 lim (A* 12 (A; x, - Sl , -s 2) -s 3 ) e -^) (17) 

A^oo 

= -si, -s 2 , -s 3 ), 

as desired. □ 

The paper is organized as follows. In Section [21 the total integrals of the purely real 
solutions of Painleve II equation are evaluated. In particular, Theorem 12.21 gives a new 
short proof of the evaluation of the constant term of the asymptotics of the GOE and 
GSE Tracy- Widom distribution functions in random matrix theory obtained in [14] . The 
total integrals of the purely imaginary solutions are computed in Section [31 In Section HJ 
we compute the asymptotic expansion of the generic purely imaginary solution up to 
0(x~ 3 / 2 ) as x — ► 00, whose total integral is studied in Theorem 13.21 In Section [5j the 
total integral of a special solution to Painleve V equation is computed, and a new simple 
proof of the constant term in the asymptotics of the gap distribution of orthogonal and 
symplectic ensembles of random matrix theory is given. Finally, in the last two sections 
the total integrals of the densities of the mKdV conservation laws evaluated for the 
Painleve functions are computed (Section EJ), and the relations to the trace formulae of 
the scattering theory for the Dirac operator are discussed (Section [7]). 



2. Purely real solutions 

A solution of Painleve II is real for all real x if and only if the monodromy data satisfy 
S3 = si, s 2 = s^. (18) 

See, for example, page 158 in [22]. The constraint on the monodromy data shows 
that if |si| = 1 then Si must be ±i and s 2 can be any real number. If |sj_| 7^ 1 then 
s 2 = (si + sl)/(l — |si| 2 ). If s 2 7^ 0, then u(x) has infinitely many poles; specifically 
([3lj; see also page 349 in [22]), for purely real solutions with s 2 7^ 0: 

fy/2 \ 
u{x) ~ i-^/xtan I — - x 3 ^ 2 + 0(lnx) ) as x — > +00. (19) 
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Since we want to integrate u(x) we will assume S2 = 0, and thus that si is purely 
imaginary. If |si| > 1 then u(x) again has infinitely many poles; specifically ([34]; see 
also page 349 in [22]), for purely real solutions with |si| > 1: 

u(x) ~ iv^x/sin f~(-x) 3/2 + 0(ln(-x))^ as x -> -00. (20) 

There are two cases of global purely real solutions: 

• The purely real Ablowitz-Segur solutions [4T1 [42] with monodromy data 

- 1 < is 1 < 1, s 3 = ~sj=-s 1 , s 2 = (21) 
and asymptotics 

3 {-xf* + (3 log(8(-*) 3 / 2 ) + + O (p^) 

(22) 

as x — > +00. (23) 



(3:=^- log(l - |si| 2 ) < 0, := -7 - argT(i/5) - arg Sl , (24) 
and Ai(x) is the standard Airy function. A representative solution with s\ = — i/2 




is shown in figure 2(a) 



The Hastings-McLeod solutions [25] with monodromy data 

s x = ±i, s 3 = s 2 = 0. (25) 
and asymptotics 



"(•r) )>•! \j — - + 0((— x) 5/2 ) as x — > -00, (26) 



2 

-(4/3)x 3 / 2 ' 
.1/4 



u(x) = isiAi(x) + O I — ) as x — > +00. (27) 



The solution with si = — i is shown in figure 2(b) 



The error estimates above come from [T4] . These solutions have no singularities for 
finite x [2j [25]. Both of these solutions look like the Airy function (up to a constant) 
as x — > +00. However, as 1 ^ —00, their asymptotic behaviors differ dramatically: 
the Ablowitz-Segur solutions decay, whereas the Hastings-McLeod solutions grow. We 
begin with the total integral for the Ablowitz-Segur solutions. 

Theorem 2.1. [Purely real Ablowitz-Segur solutions] Suppose that u(x) is a 
solution to the Painleve II equation (Qp with monodromy data — 1 < isi < 1, S3 — 
si = — si, s 2 = (that is, with asymptotics given by fiUty) and (f^j j. Then 

/;j^) d , = i.o g (I±g. (28) 
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5 10 

(a) The real Ablowitz-Segur solution with si = 




oL 

-15 



10 



-10 -5 5 

X 

(b) The Hastings-McLeod solution with si = — i. 



Figure 2. Plots of purely real solutions to Painleve II. 



Proof. Consider A approaching in the region Q 2 (see figured]). That is, set 

P(x) :=* 2 (0;a:). (29) 
Since the Ablowitz-Segur solutions are integrable on the entire real line, we could choose 



(30) 



a = — oo or a = +oo in (jlip ; we pick a = +oo. Then 
lim P(x) = lim e u (+^^C = C. 

x^+oo £— >+oo 



We therefore find C by analyzing the Riemann-Hilbert problem as x — > +oo. This 
analysis is done in [22J, Chapter 11, Section 6, so we merely provide a short sketch of 
the argument. Note that since s 2 = the jump contour consists of only four rays. We 
use the scalings 

t:=x 3/2 , $(z;t) :=V(\(z);x), 4>{z) := i~z 3 + iz. (31) 
x l l A 3 

Using standard contour deformations, this Riemann-Hilbert problem for $(z; t) may be 

transformed into the following Riemann-Hilbert problem for $> def (z;t): 

$ dcf (z; t) is analytic in C \ { 3?z = ±- 



[$ det (^;t)e 



t0(z)<7 3 



1 

Si 1 

1 s x 

1 



(32) 



2 G < Qtz 



oo. 



Here the two jump contours are oriented from — oo to +oo. Under this deformation, 



$(z;t) = $ dcf (z;t) 



def +w * M for z e o 2 n jo < < ~} . 



(33) 
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Then, from a standard Riemann-Hilbert problem small norm argument 

/ -2t/3\ 

$ dcf (z; t)e*^ )<T3 =1 + I —j=- J as t +00. (34) 
Undoing the contour deformations gives 

C= lim P(x) = lim $ dcf (0;t) ( 1 ° | = ( 1 ° . (35) 

k->+oo I Sj 1 / \ Si 1 / 

Thus 

p. . / cosh U(+oo, x) + isi sinh U(+oo, x) i sinh U(+oo, x) \ 

I — i sinh f/(+oo,x) + Si cosh f/(+oo,x) cosh [/(+00, x) / 

The analysis for x near +00 goes through even if s\ = ±ia, a > 1 (and S3 = — s\, 
S2 = 0). We will use this fact when studying the Hastings-McLeod solution below. 
However, for the analysis at x near —00 the analysis is different for s± = ±i (since u(x) 
is not integrable at that endpoint). 

The analysis of \I/i(0;x) as x — ► —00 is identical for both the Ablowitz-Segur 
solutions and the generic purely imaginary solutions. This calculation is carried out 
below as part of the proof of Theorem 13 .21 Specifically, for the Ablowitz-Segur solutions 
equation ( 1H71) holds with S3 = —si. At A = the two functions ^i(A;x) and ^(A;^) 
are related by a multiplicative jump: 



lim P(x) = lim \[> 2 (0;x)= lim *i(0;x) 



x— ►— 00 



1 

si 1 



l£3_ 



(37) 



1 



(3? 



1 -s 3 W 1 \ = / Vi - sis 3 

vi - sis 3 V _Si 1 / V Si 1 / V 

Combining (135]) and (13"7|) and using s 3 = — Si shows 

cosh Z7 (+00, —00) + isi sinh [/ (+00, —00) i sinh U(+oo, —00) 
-i sinh U(+oo, —00) + si cosh £/ (+00, —00) cosh U (+00, — 00) 

^ Tfcf 

The (21) entry gives 

(-i + s 1 )e u{+co ^ co) + (i + s^e" 17 ^ 00 '" 00 ) = 0. (39) 
Solving for U(+oo, —00) gives 

y u(y)dy = - log f ^— — J + 2i7rm (40) 

for some m G Z. Since is purely real, m = 0, which gives equation (|2"81) . □ 
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Next we compute the total integral of the Hastings-McLeod solutions. Since these 
functions are not integrable near x = — oo we will subtract off the nonintegrable part. 

The integral of the Hastings-McLeod solution with si = — i appears in the Tracy- 
Widom distribution functions that arise in random matrix theory [4"5l |4*6] . The proof of 
Theorem 12.21 below is a new, shorter way to show a result that was obtained previously 
by the first three authors using the asymptotics of orthogonal polynomials in [3] . 

Theorem 2.2. [Hastings-McLeod solutions] Suppose that u(x) is a solution to the 
Painleve II equation (QP with monodromy data s± = ±i ; S3 = =pi, S2 = (that is, with 
asymptotics given by /[2b]) and (27\)). Then, for any eel, 

jf u{y)dy+ J ^)-i Si ySj d2/ = _ isi ^ c | c |V2 + iSl I log ( 2 ). (41) 

Proof. We set s% = i. The alternate case s± = — i follows immediately from noting that 
if u(x) is a solution to (JTJ) then so is —u(x). 
Take X E Q2 and define 

P(x) :=* 2 (0;ar). (42) 

The Hastings-McLeod solution is integrable at x = +00, so set a = +00. The constant 
matrix C was computed above in (135]) in the section on Ablowitz-Segur solutions. 
Explicitly, 

C = lim P(x) = I 1 ° I , (43) 

and therefore 

( e- u ^ isinhf/(+oo,x) \ 

[ > y ie -U(+oo0) COS hf/( + 00,x) J' 1 ' 

Now we compute the asymptotics of P(x) as a; — > —00, taking into account the 
nonintegrable term using a (yf-function. This Riemann-Hilbert problem was analyzed 
in [13], and we give a sketch of the argument. Recall that the function \l/(A;x) and 
m^(X;x) used in [13] are related as in ( f!4l . Define 

5 (A) := (A 2 - I) 3 / 2 (45) 

with branch cut on [—1,1] and sheet chosen so g(X) ~ A 3 as A — >• 00. Then set (see 
(6.17) in p 



m§(A;x) := m« ^^A;xj ei^W-^-f t := ^( 



x) 3/2 . (46) 
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By standard changes of variables we can transform m 9 (A; x) to m^ 23 \\; x), which solves 
the Riemann-Hilbert problen|§| 

m^ 23 \\; x) is analytic in C \ E 

mf\\- x) = m { ™\\; x ) e -^/^ e -itg-W^ V ^A^g+W^^/\ A £ E (47) 

m (23) (A;x) = / + 0(A- 1 ), A -> oo, 

with the contour E and the constant jump V nM given in figure [31 For A G f2 2 , 

(A; x) = e imJ3/4 m (23) (A; x) e - 0(A;:c),T3 e- imJ3/4 . (48) 

In [Hj the analysis includes s 2 ^ 0, but for us S2 = S5 = I. As x — > —00 (that is, 
t — > +00), formally the jump approaches the identity on all portions of the contour in 
figure [3] with the exception of the jump S^ 1 S^ 1 = SqSi on the interval [—1, 1]. Now 




Figure 3. The deformed Riemann-Hilbert problem for the Hastings-McLeod solutions 
to Painleve II. 



nr 



'(A; x) = (I + 0(r 1/2 ))m mod (A) [H], where m mod (A) solves 



m mod (A) is analytic in C \ [-1, 1] 
m - od (A) = m° lod (A)e- i " CT3/4 (° 1 | <> 



^m mod (A) =/ + 0(A~ 1 ), A ^00. 
This problem is solved explicitly by 



i7TfT3/4 



AG [-1,1] (49) 



m mod (A) 



/(A) := 



-17TCT3 



/4 { /(A) + /(A)- 1 /(A) - /(A)- 1 \ _ 3/4 



/(A) - /(A)" 1 /(A) + /(A)- 1 



A- 1 
A + l 



1/4 



(50) 



We use the notation m^ 23 ' (A; x) to correspond to reference [14] 
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Here /(A) has its branch cut on [—1, 1] and /(A) ~ 1 as A — > oo. Undoing 
the transformations from \J/ 2 (A;x) and using /+(0) = e 171 "/ 4 , g+(0) = — i, and 
lim m (23) (A; x) = m mo (A) outside of small neighborhoods of ±1 we have 

x^— oo 

lim P(x)e ta3 = lim * 2 (0; x)e i * 9+(0)<T3 = e i7rCT3/ V™ od (0)e- i7RT3/4 

x— ►— OO X— ►— oo 

= v^i ij ( 51 ) 

Combining (jH]) and (l5Tj) gives 

/ e -c/(+oo )a: )+t i s i n hf/(+oo,x)e- i \ _ ( 1 i . . ro 
x ™oo 1 ie -c/(+oo,x)+t cos hf/(+oo,a;)e^ J ~ T" I i 1 1 ' 



The (11) entry is equivalent to 



y w(y)dy + — (-x) 3/2 j = — . 



Hence, for any fixed c, 



tr + OO PC ^2 \ 

y u(y)dy + J u(y)dy+ — (-x) 3/2 J = — . 

From the asymptotics of tt(y) we see that u{y) + ^J\y\/2 is integrable at — oo, so 



(53) 



(54) 



X— > — oo 




lim exp I / u(y)dy+ j \ u(y) + \/ — 1 dy - — c\c\ 1 | = — . (55) 



Taking a logarithm shows 



-OO PC 



c J — oo 



w(2/)dy + / I u (y) - y ) dy 




y/2 1 
— isi— c\c\ l/2 + isi- log(2) + 2i7rm 



(56) 



for some m G Z. Since is purely real, we see m = 0, which shows (|41l) . □ 
3. Purely imaginary solutions 

Solutions to ([1]) are purely imaginary if and only if (page 159 of [22]) the monodromy 
data satisfy 

s 3 = -si, s 2 = -sjj. (57) 

All purely imaginary solutions are global (page 297 in [22J). There are two distinct 
asymptotic behaviors: 
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• The purely imaginary Ablowitz-Segur solutions [HJ H2] with monodromy data 

si ef, s 3 = s 2 = (58) 

and asymptotics 

//(,■) sin ( -(-xfl 2 + -d 2 log(-x) +A+0 ^ log ( - ^ 



x)V4 V3 4 bV 7 7 V(-^) 5/4 / (59) 



(61) 



as x -> — cx), 

/ e -(4/3)*»/a 

m(x) = isiAi(x) + O I ^ — ) as x — » +oo, (60) 

with 

d 2 := -log(l + | S! | 2 ) , d>0, 

7T 

3 l2l , n , 7T ^ f d 2 \ 

(f) := -d log(2) - - - argT I i— J - args x . 

The generic purely imaginary solutions [27J [T3] with monodromy data 

9f(si)^0, s 8 = -5T, ^2 = ^^ (62) 

and asymptotics 

"(■'•) = ; ^4 „ sin I -(-x) 3/2 + -d 2 log(-x) + 0| + ' 1<,R( 



x)V4 V3 4 bV 7 V V(-^) 5/4 /' (63) 

x — ► — oo, 

w(x) = ier\[^-+ ; cos ( ^^x 3/2 - -p 2 logx + 6M + O , 
1 ; V 2 (2x) 1 /4 ^ 3 2^ y W (64) 

x — > +oo. 

Here d and are given by (IBTj) . and 

P 2 : = -- log(|s 2 |) = - log ^,tJ Sl || . P > °> 

a := -sgn(S(s 1 )), (65) 

37T 7 

:= -— - -p 2 log 2 + ar g r(ip 2 ) + arg(l + s 2 ). 

Note that the asymptotics as x — > — oo are exactly the same for both types of purely 
imaginary solutions. We first find the integral of the purely imaginary Ablowitz-Segur 
solutions. The result is the same as for the purely real Ablowitz-Segur solutions. 

Theorem 3.1. [Purely imaginary Ablowitz-Segur solutions] Suppose that u(x) 
is a solution to the Painleve II equation (TJP with monodromy data Si 6 I, S3 = — Si, 
s 2 = (that is, with asymptotics given by / f5Pj) and / fgQj) ). Then 



/ +oc 1 /l + isi \ 

u(y)dy = - log f - _ — j = iarctan(si) = iarg(l + isi). 



(66) 
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Proof. The asymptotic analysis of the purely imaginary Ablowitz-Segur solutions is 
exactly the same as that for the purely real Ablowitz-Segur solutions. The proof of 
Theorem 12.11 applies without change through equation ( HO"]) : 



/+°° i / 1 _|_ j s \ 



(67) 



Assume S3 = — si and S2 = and parameterize the purely imaginary Ablowitz-Segur 
solutions si) by Si. Note that si = corresponds to to the solution u(x; s\ — 0) = 0, 
and in this case clearly m = 0. We now show continuity of the total integral u(x)dx 
with respect to Si for Si G JR., which shows m = in fIBTj) . The Fredholm theory 
for Riemann-Hilbert problems shows that, for fixed x, the solution \I/ to ([3]) is either 
meromorphic in s% or there is no solution for any Si ([22] Corollary 3.1). Furthermore, 
the associated Riemann-Hilbert problem has a global solution for all Si G M assuming 
S3 = — Si and S2 = (see [22] Theorem 5.6 and note the condition in (5.5.1) should read 
l s i + S3 1 < 2). Combining these two facts shows that \1/ is analytic in si for the purely 
imaginary Ablowitz-Segur solutions, and thus u(x; si) is continuous in si. To show the 
total integral is continuous in si we fix L > large and show 

lim ( / (u(x; Si) — u(x; s[))dx + / (u(x; Si) — u(x; s[))dx 

1-iW-oo (6g) 

(w(x; Si) — u(x; s[))dx J = 0. 

The continuity of u(x; Si) with respect to Si shows the limit of the second integral is 
zero since the region of integration is compact. For the third integral, use (1601) to write 

u(x; Si) = isiAi(x) + E + (x; Si), (69) 

where \E + (x; Si)| < B + (x) for some B + (x) G L 1 uniformly for Si G [s^ — e, s^ + e]. So 

lim / (u(x; si) — w(x; s' 1 ))dx = lim i(si — s^) / Ai(x)da; 



si 



H-OC 

lim / (£+(x;si) -E + (x;s' 1 ))dx (70) 



= 

by the dominated convergence theorem. For the first integral, use fl59l) to write 

u ^ = (_^/4 sin Q(~ x)3/2 + i rf2 io s( _x ) + ^) + si )> 

where \E~(x; S\)\ < B~(x) for some B~(x) G L 1 uniformly for S\ G [s^ — e, s\ + e]. 
Direct computation shows that 

£ sta (I<-) 3/2 + ^ M-*> + *) d * 



^ -i0r3(l-id 2 )/4 

3 
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where E n (z) := ^-dt has a branch cut in z on (— oo,0). The right-hand side of 
fl72|) is continuous in si. Therefore, by the dominated convergence theorem, 

lim / (u(x; Si) — u(x; s[))dx = 0. (73) 

This verifies equation ( 1681) . □ 

We now compute the integral of the generic purely imaginary solutions. The 0{x 1 ^ 2 ) 
term in the asymptotic expansion (1641) as x — > +oo is not integrable, so we will subtract 
it off as in the Hastings-McLeod case. The 0(x -1 / 4 ) term is integrable because of the 
cosine factor. However, the 0(x _1 ) term is not integrable, so it must be computed and 
subtracted off as well. The explicit form of the 0{x~ v ) correction to the asymptotics 
( |64l) was formally calculated via the analysis of a certain nonlinear integral equation 
equivalent to ([TJ) in [32]. The asymptotic expansion for u(x) up to the 0(x~ l ) terms 
turns out to be 

/ \ • fx , i^P ( %V2 3/2 3 2 \ 3iap 2 

uix) = ia y2 + j2xW cos [— x ~2 P logx+e )- — 

(74) 

■cos ( 2 x 3/2 - -p 2 logx + # 



4a; \ 3 2 



+ 0(x 3/2 ), x -> +00. 



With the first two terms already known, the third and the forth terms of this formula 
(and, in principal, the terms of an arbitrary higher order) can be formally derived via 
substitution into the Painleve equation (0Q) (or to the nonlinear integral equation of |32j). 
It should be emphasized that even the formal derivation of (|74|) is quite challenging; 
indeed, because of the presence of the growing term yVx/2, it is much more difficult 
than the similar derivation of the correction terms to the semi-linear asymptotics (j22p . 
A serious additional question is the justification of the asymptotics (j71|) which can be in 
principal done using a priori information of the structure of the asymptotic series which 
in turn can be extracted from the Riemann-Hilbert analysis (compare to the approach 
of [15]). In Section H] we will present an alternative and rigorous derivation of (IT4"|) using 
the direct asymptotic analysis of the Riemann-Hilbert problem ([3]). It also should be 
noticed that, in fact, we do not need to know the 0{x~ x ) terms a priori in the proof 
of Theorem 13 .21 The 0{x~ x ) term that must be subtracted off to make the integral 
finite arises naturally during the computation. However, note that the oscillatory term 
of 0(x _1 ) in (|74l) will not be subtracted off because it is integrable. 

Theorem 3.2. [Generic purely imaginary solutions] Suppose that u(x) is a 
solution to the Painleve II equation with monodromy data $s{s\) 7^ ; S3 = —~s±, 
$2 = (si — + |si| 2 ) (that is, with asymptotics given by (Ofjj) and (64\)). Define p 2 
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as in (E3Jj. Then, for any c > 0, there exists m G Z such that 



u(y)dy + I ( u(y) - io-y| + i<7^- ) dy 



ia{ arg(l + iasa) - ^- log 2 + arg ^ Q + i^ ^ 



(75) 



3/2 3p 2 , 

H c 3/i ^-logc + 2vrm 

3 4 



where a := —sgn(Q(si)) and T denotes the Gamma function. 

This result determines the total integral up to an additive factor of 27rim for some 
m G Z. 

Proof. Since the solutions are integrable for x near — oo, pick a = — oo and consider 
U(—oo,x). For convenience we consider ^i. Set 

P{x) :=*i(0;x). (76) 

Now we compute 

C= lim P(x) (77) 

X— > — oo 

using the methods in [22]- Start with the scalings 

A 



z : 



t := f-x) 3 / 2 , 



l"^) 1/2 ' ^ ( 78 ) 

t) := *(A(z); x), 0(z) := i-2; 3 - Lz. 

The solution can be transformed using standard algebraic manipulations to 

ty def (z;t) which solves the Riemann-Hilbert problem on the deformed contour shown 
in figure HI wherein 

*-(*:)• X)- 

The normalization for the deformed Riemann-Hilbert problem is 

# def (z; t)e^ z)ta3 =1 + Oiz- 1 ) as z -> 00. (80) 

In particular, 

^(A;x) = ^ def (^;t) i-'i-a i" 51 ^ 
\ -si 1 

for A G fl {3A > 0} fl {|A| < 1/4} when x < —1. The jump matrices for 
\l/ def (z; t)e t ^ (T3 off the real interval [— |, |] decay to the identity as t — ► 00. Indeed, 
from page 328 in [22], 

$«faf = $i'(/ + 0(r 1 / a )), (82) 
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\S2Sn 1 



So 
— s — 




Figure 4. The deformed Riemann-Hilbert problem for the generic purely imaginary 
solutions and the Ablowitz-Segur solutions to Painleve II as x — > —00. 



where ^> D (z]t) solves the model problem 



^ D (z;t) is analytic off 



1 1 

2' 2 



^ D (z;t) does not have a non-square integrable 
singularity at the endpoints ± - 



*°(z;t)=*°(z;t)S D , ze 
$ D (z;t)e te ^ ->/, 



1 1 

2' 2 



.4 



This problem is solved by 



V D (z;t) 



m 



/(*) 



z — > 00, := i^-z 3 — iz. 



z + 



1 \ M 



2 



-1 



27Ti 



(83) 



The function f(z) is defined with its branch cut on [— |, |] and satisfies /(z) 
z — > 00. It follows that 



(84) 

(85) 
1 as 



/ 1 S3 

lim P(x) = lim #i(0;x) = lim ^J f (0;t) 1 " SlS3 1 " SlS3 



-S3 



(86) 



-si 



Using / + (0) = VI - S1S3 gives 

C = lim *i(0;x) = 



1 



-S3 



(87) 
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We note specifically that this gives 

(P(x)) n + i(P(x)) 21 = iZ^- e g (-^. (88) 

V 1 + \ s i\ 

Now we analyze P(x) as a; — > +00. This limit does not exist since is not 
integrable at +00. However, the limit of P(x) times an appropriate decaying factor will 
exist. In [22J (see page 346) it is shown that 

(A; x) = (I + 0(x- 3/4 ))$(z; x), (89) 

where ^(z; x) is the solution to a model Riemann-Hilbert problem. Let be \I> in the 
region flj.. Then, by the computations in 



*l(0;x) = ( ^ ^ j e -™/4 2 7^ 3 /4 (2 g ) -<x3/2 



X 



where 



z '■— ~^u2i t := x 3 / 2 , v:=- log ias 2 , Q := iT (1/ + 1 1 2 . 

X 1 ' 2 17T a/27TS 2 

1/i " (91) 



C(*) := 2 A /it^ - itsgn(^)^ ^ + 1 

and Z RH (C) is a function built out of parabolic cylinder functions. Specifically, for 

Z ( ° 2 ^ i^.-iOC) ±D„{Q ){ 1 ) Q ' (92) 
where D u (() is Whittaker's parabolic cylinder function satisfying 

' ft + (,4-?)ft = l (93) 



dc 2 V 2 4 

Note that \I/i(0; x) is uniformly bounded independent of x. From Whittaker and Watson 
ESJ, Section 16.5, 



f 1 Z 7 12 ) := lim Z™(C( Z )) 

/ r(|) 2 - y /a-ya r (|) 2 "/2 \ / i7r(l/+1)/2 \ 
_ 9 -a 3 /2 r(i+|)^ r(p|j^ / e ^ + )/ \ /2 

Equation ( 1901) implies 

(*!(0; x)) n + i($!(0; x)) 21 = ^l^i 1 + °) Z n + 5(1 - <r)Z 21 ], (95) 
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where 

5 := e -^/ i 2 7v l\2Q)- 1 l 2 e^ tV ~ 2l h vl2 . (96) 

Assume for the moment that a = +1. Then 

($i(0; x)) n + i($!(0; x)) 21 = 2ie inu/i 2~ 7u ^Q 1 / 2 Z 11 e itV2 / 3 r u/2 . (97) 

Now from (1591 and using the fact that \l/i(0; x) is uniformly bounded independent of x 
we find 

lim (P(x) n + iP(x) 21 )e- ix3/2V2/3 x 3u/A 

= lim ((^ 1 (0;x))n + i^i(0;x)) 21 )e~ ix3/2 ^ 3 x 3 ^ 4 (98) 

= 2ie i ™ /4 2- 7l/ / 4 Q 1/2 Zn- 
Along with ( 1551) this gives 

lim exp (u(-oo, x) - ix 3/2 V2/3 + (3i//4) log(z)) 

' 2ie i^/4 2 -7,/4gl/2 Zi ^ 



1 + Isil 



(99) 



1 — isi 

Writing U(—oo,x) = £/(— oo,c) + U(c,x), (1551) implies that 

aS+oo^f/ u (y) d y + J c + i^) dy ~ + T logc 



1 + Isi 



12 



(100) 



; 2ie iW4 2 - 7l V4gl/2 z 



1 — isx 

for any c > when a = +1. It follows that the integral (u(y) — iy/y/2 + 3u/Ayj dy 



is convergent, and hence 

I n,(iA — \* I 

2 Ay 



exp ( / u(y)dy+ I [ u(y) - uj | + ) dy 



y/1 + | gl | 2 (1 + gigg) T(l + Z/) 3W4 9iV4 W2 ^ c 3/2 _ ^ 

i-i Sl 2V2 S2 r(i + |) e e exp ^ 3 c 4 iogc y 

1 + T(l + u) 9u/4 f.V2 3/2 3i/ ' 

—2 ' exp i— c ' -— logc 



(101) 



2V2(1 + | Sl |2)l/4| Sl _ T(l + |) M 3 4 



- (l + i^rQ + f) 2 -W4 exp ^ c s/2 _ 3, los ; 

" (27r)V2(i + | Sl |2)i/4| Sl _sl|i/4^ ex P ^ 3 c 4 1o § c 

The first equality follows from the definitions of Q and Z\\ in (1911) and (1941) . respectively, 
the second follows from s 2 = [s\ —~s\)l (1 + l s i| 2 )> an d the third follows from the identity 
(see (6.1.18) in [3]) 

r(2*) = -^2 2 T(z)r + ~Y (102) 
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The fact that the right-hand side of fllOip has modulus 1 follows automatically from 
the fact that u(x) is purely imaginary. However, this can also be checked directly using 
(6.1.29-31) in [3]. Equation ( 1751) with o = +1 follows by taking the logarithm of both 
sides of fllOip and setting v = — ip 2 . 

Now assume a = — 1. This result can be obtained from the a = +1 case via Lemma 
11.11 We also give a direct proof as follows. Using the definition of Z 2 \ in (1941 . 

($!(0; x)) n + i($ 1 (0; x)) 21 = 3pL 2-VV^ +2 )/ 4 2 5 ^V if ^ 3 r/ 2 . (103) 
From (1591 and the fact that ^1.(0; x) is uniformly bounded independent of x, we see 



lim (P(x) u + iP(x) 21 )e lx / V2/3 x 



ia; 3 / 2 v / 2/3^-3i//4 

lim ((^^x^n + K^i(0;x)) 21 )e ix:i/2V ~ 2 / 3 x-^ 4 

->+oo 

£^l]_2-l/2j e i7r(i/+2)/425^/4 



(104) 



1 V2 ' 2 



From 



lim exp (u(-oo, x) + ix 3/2 V2/3 - (3z//4) log(x)) 



1,112 p/l'N (105) 
I 1 1 L \ 2]_2~l/2- i7r(y+2)/425^/4 



r(| + | 



Writing Z7(— oo,s) = [/(— oo,c) + [/(c, x), (11051) shows that 



i + h 2 r(-i 



2- 



. 2 -l/2 ie i7r(v+2)/4 2 5i./4 (1Q6) 



r(|- 

= (2ny/ 2 (l + \ Sl \ 2 )^\ Sl -^ W4 

(i - iSl )r(| + 1) 

for any c > if cr = —1. Therefore the integral ^u(y) + i^y/2 — 3z//4yj dy is 



convergent, and hence ( 1751 with o = — 1 follows by taking logarithms and using 
v = ip 2 . □ 

4. Direct computation of asymptotics of u(x) in the generic purely 
imaginary solutions 

In both [14] and [22] the authors write down asymptotic expansions of the purely 
imaginary solutions to the Painleve II equation for large positive x. In this section 
we calculate the higher order terms for these expansions. Specifically, we will calculate 
the 0(x _1 ) terms in the asymptotic expansion (16~41 of the generic purely imaginary 
solution as x — > +oo and show: 
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Theorem 4.1. Let u(x) be a generic purely imaginary solution of the Painleve II 
equation (QJ) with asymptotic expansion $U3j) as x —* — oo. Then 



u{x) 



,at/ 2 + (2a:)V4 



lap / *- v * 3/2 
cos — — X ' 



'2V2 



+^-cos ( 2 



4x 



3 

2^2 



. 3iap 2 
2 -p log* + - — 



., x 3/2 - -p 2 log 2 + 6> 



(107) 



-3/2\ 



-OO, 



where a, p, and 6 are defined in /[U5\) . 

Much of the notation is inherited from 



. We note that 



v = lp 



t 



X 



3/2 



Here the solution u(x) to (IT]) is obtained as 



u[x) 



i(T \/? + 2 ^ lim i z Xn{z)) 
y z z^oo 



(108) 



(109) 



where Xui 2 ) is the 12 entry of the 2x2 matrix valued function that solves the ratio 
Riemann-Hilbert problem 

X is analytic in C\7 

x{z) — > I as z —> 00 (HO) 
X+(z) = x_(^)G(2;) on the contours 7. 




Figure 5. The contour 7 for the Riemann-Hilbert problem for x( z )- 



The jump G(z) is given in [22] (9.5.61) and (9.5.62). As illustrated in figure [5] the 
contour 7 is the union of the several contours ji, i = 1---6 and C m , m = 0,u,d. 
The contours ji are the anti-stokes lines and the contours C m , m = 0,u,l, are small 
circles oriented clockwise around the origin, +i/v^2, and — i/v2 respectively. On each 
of these contours the jump G(z) has a different definition and we write G(z) = Gi(z) 
or G(z) = G m (z) to denote the corresponding jump on each contour. As in (9.5.73) in 
[22] one can write 

lim (2x12(2)) 

= -^[f (°( z ) - 7 )i2^ - J (( X -(z)-I)(G(z)-I)) 12 dz. 
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To obtain (I6"4"|) the authors of [22J proved that equation fill II) reduces to 

lim (zxn(z)) =~ f G (z) 12 dz + 0(x~ 3 / 2 ). (112) 

See (9.5.74) in [22J and note "Stv = 0. It is exactly the 0(x -3 / 2 ) terms that we wish to 
now compute. To calculate these terms there are several things to check. The following 
three assertions, once proven, will establish the desired result. 

(i) In [22J , the authors do not compute / G (z)i 2 dz explicitly, but they compute 

2m Jc 

1 f - 

the integral / Go(z) 12 dz where Go is an approximation of Gq. The error 

2m Jc 

from using this approximation is written as 0(x~ 3 / 2 ), which could contribute to the 

0{x~ 1 ) term in u(x) (note the \fx in (11091) ). However, the error in the off-diagonal 

entries is actually higher order and does not contribute to the 0(x~ l ) term. 

1 f 

(ii) The contribution to u(x) from the integral / ((x-{ z ) ~ -0 (G(z) — I)) 12 dz is 

27ri y 7 

+ % cos ( 2 M MO - 4^/3 - 26). 

(iii) The contribution to u(x) from / G(;?)i2dz is broken down into the sum of 



the integrals on each component of the contour 7. In [22] it is shown that 
||Gj(z) — I\\ < cexp(— (2rc) 3 / 2 1 ^| 2 ) for some constant c. Consequently, the integrals 
j ' (Gi(z) — I)i 2 dz will not contribute to the term we wish to compute. The 
contribution from the integral on Cq is handled in the first assertation. The 

contribution from / G u (z)i 2 dz and / Gi{z)v>dz is — rr. 

2niJ Cu uK m 2niJ Ci lK m 2 * 

To prove assertion 1 we first proceed to verify the following Lemma: 
Lemma 4.2. For z G Co, the jump Gq has the expansion 

Go = Bo(z)e-*F°* (Mo + Mf) e^(B Q (z)y l + O (r 3 / 2 ) , (113) 

where 

f )■ ^==(T -40- (114) 

The function B (z) is analytic in a neighborhood of the origin and is defined as (see 
JEM/) 

{Y{z){C{z)8{z)r R(*)>0 
B (z) = < „ (115) 
[Y(z)iaa 1 (ias 2 )-^{C{z)6(z)Y 3 $t{z) < 0. 

27ie functions Y(z) and 5(z) are given by 

~ 2 [ a{p - /T 1 ) + ' 



5 2 (^) = I (^±i^_^ ) ( V= L i n(iaS2 ) 



(116) 



(- 2 + l) 5 + *. 



17T 
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and ((z) is as in equation /jyi\). 

Proof. Consider (9.5.50) and (9.5.43) of [22]. Combining these two facts gives: 



* = Bo(z)e-*P°* \\\\Q- 



2 



From (9.5.55) and (9.5.56) in [22] we can then obtain 



Q^eV 4 



(117) 



m„J i + ^ + o(r 4 ) ^ + o(r 3 ) 



Go-^-We-'. "( Q + g ^ + 0(c _ 5) 1 _^ + 0(c _ 4) y (iig) 



where 



C 

From the definition of ( 3 in (19 ip one can see that C 3 = ct~2 z 3 (1 + 0{z)) for some constant 
c. Equation CI 1 1 3 [) then follows. □ 

Consequently, we have that 

G (z)dz = ! G (z)dz + ! G { o\z)dz + O(t-^), (119) 

Co J Co J Co 

Go := B o (z)e-^M o e^(B (z))-\ 
£$> := B (z)e-^M^e^(B (z)r\ 

The first integral on the right hand side of equation fl 1 1 9 [) is computed explicitly in [22] 
and is what gives rise to the cosine term of the expansion ( 164"1) . The second integral of 
(11191) does not contribute to the next order term of fl 1 1 2 j) . moreover: 

(Gf) 12 (^)dz = 0. (121) 

Co 

To show this, we first write down the (12) entry of Gq using the definition of Mq . 
Simple algebra yields that 



(G i * ) ) 12 (z) = -~— ^ -f \ R(*)>0, 



2 C 2 



r - (122) 

2 

The function /? has a branch on the imaginary axis. However, since B is analytic in a 
vicinity of the origin we can deform the contour of integration through the branch so that 
it does not pass through the interior of Co- However, as a function of z, ( 2 is analytic 
in a vicinity of the origin, but has a zero of multiplicity two at z — 0. Consequently, 
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when we deform the integral we pick up a residue from the origin. From the definition 
of ( in (JHD we can write C 2 = -M^z 2 (l + 0{z 2 )). With this it is clear that 



a -it2V2.l Co z> tj2iz v "»=» 

The last inequality is due to the fact that /3(0) + /3~ 3 (0) = 0. This proves the first 
assertion. 

Next we check assertion 2 by showing 

1 



2?ri 



(123) 



((*-(*) (G(z)-I)) l2 dz 

= -^7 + ^7 cos f 2 M M*) - 4 ^/3 - 2fl) + O (r 3/2 ) . 
ot ot V / 

Let C be the Cauchy operator on L 2 ^) defined for / £ £ 2 (7) by 

(C7/)(z) := / -^ds for * $ 7. (124) 

27T1 J 7 S - Z 

Let C_ denote the boundary limit of C defined for / £ L 2 {j) by 

(CL/)(s):= lim(C7)(A (125) 
where z' is on the right-hand side of 7 and z £ 7. Define, for / £ £ 2 ( 7 ), 

(C G ^jf)(z):=C4(G-I)f](z) forz£ 7 - (126) 

Note that WCg-iWl^^^l 2 ^) < I|C-Hl 2 ( 7 )^l 2 (7)||G' — ^||l°°(7)- Using the fact ((9.5.69) 
in [22]) that the jump matrix G satisfies 

Q 

I \G - J| | L 2( 7 ) ni oo( 7) < — (127) 
for some constant c we have 

||Cg-/||l 2 (7)^l 2 ( 7 ) < —q (128) 

as C_ is a bounded operator on L 2 {j). 

Suppose /i(s) satisfies fj, — I £ £ 2 ( 7 ) and (1 — Cg-i)(h — I) = Cq-iI- Note that 
from (11281) . 1 — Cq-i is invertible for all x sufficiently large. So, for x sufficiently large, 

— -^IU 2 (7) < 11(1 — Cg-/) _1 Cg-/^||l 2 (7) < 2||C G _//|| L 2( 7 ) 

<?\\r 11 nr m < c (129) 

S ^H l '-||l 2 (7)^l 2 (7)|| ( - t — 1 I|l°°(7) S — 

From standard Riemann-Hilbert theory (see, for instance, pj), for z ^ 7, 

, , r 1 f u(s)(G(s) - 1) . T 1 [G(s)-T , , , 



where 



BW ;= w (, M -/KGM-/) ds = _ /)(G _ 7)] w (131) 
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Defining E-(z) := C[(ji - I)(G - !)](*), we have 



J E-{z){G{z) - I)dz 



< c||^-||i2( 7 )||G — /||l 2 ( 7 ) 



< ||C_|| L 2( 7 )_ +L 2( 7 )||/i — J|| L 2( 7 )||G' — J|| i0 o( 7 )||G — J|| L 2( 7) 

c 

< — T". 

- t 3 / 2 



(132) 



Therefore 
1 

_ 2vri 



(X-(z)-I) {G{z)-I)dz 

i r /" HfWK^W) dsd , + (r3/2) (133) 



(27ri) 2 J 7 J 7 s — z 

As (see (9.5.64) and (9.5.65) in [22]) 

IIC - ^IU 2 (7\C )nL°°( 7 \Co) < 7 (134) 

for some constant c, in (11331) we can restrict the integrals to Co without adding a larger 
error term: 
1 



27Ti 



{X-{z)-I) (G(z)-I)dz 

1 f f (G (s) - I)(G (z) - I) 



(2vri) 



J Co J C S- Z- 



(135) 



(137) 



where Go denotes G on Co. Furthermore, as (see (9.5.66) in [22J) 

m c 
\\G - G \\L 2 (Co)nL°°(c ) < ^, (136) 

where G is defined in (9.5.58) of (22], we can also replace G by G : 
~ j (X-(z)-I) (G(z)-I)dz 

= 1 / / (Go( g )-/)(G (,)-/) 

Now we evaluate this double integral explicitly. The function Gq(z) is defined ((9.5.58) 
of 1221) by 

/ 1 v c -i2V2t/3 \ 

G (z) = B (z) . 2V - «W B^*), (138) 

V cw e 1 / 

where Q and v are constants given in (1911 above and and -Bo(^) (see (11151) ) are 
holomorphic in Co. The error term O(z) in (9.5.47) of [22] is actually 0(z 2 ), and the 
function £(z) satisfies 

= e - i7r / 4 v / t2 5/ S(l + 0(^ 2 )), z ~ 0. (139) 

The function -Bo(^) satisfies ((9.5.53) in [22] ) 

S„(0) = -L [ / ] e -™/4 r <73/2 2 7,<7 3 /4 ( 140 ) 

y2 \ ic 1 
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where a = — sgn(Q'Si). Note that G (s) has a pole at s = 0, with the residue (see 
(9.5.76) in [22]) 

A : = ResG(^) 



2=0 




n " p -2h/2t/3 \ 

o ^° (0) (141) 

ia(qt~ u e 2y ^ lt / 3 — pt u e~ 2y ^ u / 3 ) qt~ u e 2y ^ u / 3 ^1^-2^/3 

g t -u e 2V2it/3 _|_ p t u e -2V2it/3 _j (J ^^-^ e 2v / 2it/3 _ p^ e ~2 v / 2it/3^ 



o-3/4 /= mu/2 

P = ^ L, r( q = p. (142) 

(1 + SiS2)T(l/) 

There are a few typographical errors in (9.5.76) and (9.5.78) of [22]. In (9.5.76), the 
diagonal entries of the middle matrix which is conjugated by -Bo(O) should both be 0. 
The diagonal entries of the middle matrix in the last equality of (9.5.76) should also 
both be 0. In (9.5.78), T(u) should be replaced by T(—u). 
From the residues at s = z and s — 0, 

1 f f (G (s) - I)(G (z) -I) dsdz 



^JcJc s-z_ (m3) 

~( (G ( Z )-ifa z + ^ [ WM-V dz+otr*). 

2tti J c 2vri J Co z 

Noting that the (12) entry of 

(Go (z) - J) 2 = -f-J (144) 

is zero, we see the second integral on the right-hand side of (11431) does not contribute 
to the evaluation of u(x). Next, another residue calculation shows 

j_ r A(G (z)-i) dz 

2vri J Co z 

- 1 ! l AB(z)( ° S^Ms^d, (145) 
" 27i J Co V 2ABq[z) \ f e 2i ^/3 o ) B0{Z) dZ 

= A[B> (0)B (0y\A}, 

where [M, N] := MN — NM. From (j!15p . a direct calculation shows that 




(146) 
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.4 



la 



1 



A 



or* 



12 



>/2 \ 1 

^-iz/2~ 3/2 + g 2 e 4V2it/3-2vlog(i) + p 2 e -4V2it/3+2ulog(t) 



(147) 



+0 (t z 

~ + ^- cos f2|z/| log(t) - 4V2t/3 -26) +0 (V 
St 8t V /V 



where 6 1 is defined in ( 1631) . The last equality in (j!47p uses 



r(i + i y ) 



ny 



for |/ > and S2 



si - Si 



(148) 



sinh?/ ' 1 + 

We proceed to check assertion 3. The jump matrices G u and G; are written as 
\E fUl I f ~ 1 and ^l/'^ -1 . G u and G/ satisfy the symmetry 



Oi 



G t (z)dz = -a 2 ( / G u (z)dz a 2 . 



C u 



Using this symmetry we note that 

/ G u {z) 12 dz + Gi(z) 12 dz = G u (z) l2 + G u {z) 21 dz. 

J Cti. J Ci J Cn. 



(149) 



(150) 



This allows us to work only with G u . 4f and are explicit functions (see (9.5.24) and 
(9.5.34) of [22]) the later being constructed using Airy functions (see (9.5.30) in [22]). 
Using the leading and second order asymptotics of the Airy function and its derivative 
(see Abramowitz and Stegun [3]) one can write: 



G u 



/+ Y6°*{is 2 

4(2 



(0-3/2 



Ci + d\ —c\ + d\ 
ci - di -(ci + di) 



(151) 



x(/+o(icr 3 )), 

where Y(z) and 8(z) are defined in (11161) and 



2 2 7T1 / o 

£ = 2st3e"~ [ z + - 



Cl 72' rfl 72' 



(152) 



The leading-order contribution from the integral j c (G u (z) — I)dz will come from 
integrating 



-W 3 (is 2 ) 



^3/2 



ci + d\ — ci + c?i 
ci - c?i — (ci + d\) 



is 2 



^— (7 3 /2 J— <J3y— 1 



(153) 



In particular we are interested in the sum of the (12) and the (21) entries of this matrix, 
which can be written out as 

i 



(12) + (21) 



16xt(z2 + i/2)l 



(6 2 (is 2 ) - r 2 (i S2 )- x ) . 



(154) 
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We expand S(z) for z near i/\/2, giving 
5 2 {z) = — { 1 + 2ie 7ri / 4 z/(z - i/^2)3 + 8z/V i/2 (^ - i/v^) 



+0(^-i/V2)i 



iV2 | - ,2V2V 2 -7fJ +0 U 2 "^ 1 Ul57) 



Inserting the expansions (j!56p and (11571) we obtain 
(i S2 )5(z) 2 - (is 2 )- l 5- 2 (z) 



(z* + 1) 



n 3/2 



-dz 



, (2x e f^-^)Uo((z-^)§ 



3 



3 / _ i \ „// i 



x I 1 -= [z- — )+0\[z--=\ \\dz 



2 N 



n il si r 1 
cri 4 e 4 z/ / 1 



i^2> Jc u z-^ 



-dz 



(155) 



Consequently, 

(is 2 )<5 2 - (isa)" 1 *" 2 = a (2^e^u(z - i/y/2)* + O ((z - iv 7 ^)) . (156) 

Additionally, 



(158) 



= — 8crz/7r. 

Using the definition of G u , (|158p . and (11541) we have that: 

-^-[ G u (z) 12 dz - -L / G^)i 2 d^ = -^ + 0(a;- 3 ). (159) 

27T1 y Cu 27T1 J Cl 4, X 2 

Together, assertions 1, 2, and 3 and (164")) establish Theorem 14.11 
5. The GOE and GSE sine-kernel constants 

Define J to be the interval (—1, 1). Let K' x ' be the integral operator on L 2 (J, dz) with 
kernel 

K^(z,z';x) := s ™f z -*') _ (160 ) 



smx{z — z 

71 (z — Z 1 



Also let K± be the integral operators on L ((0, 1), dz) with kernels 
„(a.) 1 /sin x(z — z') , sinxfz + z') 



± vr V ^ - z + z' y v; 

Then 

P (x) := det(l - K (x) ) (162) 
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is the limit (as N — > oo) of the probability that an N x iV matrix drawn from 
the Gaussian Unitary Ensemble has no eigenvalues in f — after proper scaling so 
that the mean spacing of eigenvalues in the bulk is normalized to 1. Also define the 
determinants 

D±(x) := det(l -Ki x) ). (163) 

Then D + (x) and \{D + {2x) + D_{2x)) are, respectively, the limits of the probabilities 
that a matrix drawn from the Gaussian Orthogonal or Gaussian Simplectic Ensembles 
has no eigenvalues in (0, ~) after scaling so the bulk spacing of eigenvalues is normalized 
to 1. Dyson [16] conjectured and Ehrhardt [18] recently proved that 

Theorem 5.1. 

^ x 2 x logx log 2 , log 2 3 ... „. ,,„,, 

\ogD ± = =f —+-2—±-2—+-C (-l)+o(i) as x^ +oo, (164) 

& 4 2 8 24 4 2 v ' K ' 

where ( is the Riemann zeta function. 

The o(l) terms are given by an explicit, asymptotic series. We give a short alternative 
proof of this theorem that will follow from Lemmas 15.21 and 15.31 To begin, we express 
logD±(x) in terms of the definite integral of a solution to the Painleve V equation. This 
function arises in the solution of a Riemann-Hilbert problem (studied in [12]) that is 
associated with the sine kernel. Set J := (—1, 1) and let m(z; x) satisfy the Riemann- 
Hilbert problem (see (1.11) in [T2j) 

' m(z; x) is analytic for z J 

g2ixz 

m + (z;x) = m^(z;x) ( _ 2ixz o ) on J 



2 r (165) 



m(*;*) = I + (!]as™ 



z t 

Here J is oriented left to right. Define m\(x) by 

m{z; x) = I + + O (\) a s z -> oo. (166) 

z \ z J 

Then set (see (4.31) in |T2] ) 

i{x) := 2i(mi(x)) 21 = -2i(m!(x))i 2 . (167) 

It is shown in [12] that £(x) is related to a solution of the Painleve V equation. Indeed, 
let u(x) be the regular at x = solution of the Painleve V equation 

d 2 u _ f du Y 3m -1 + 2u{u + l) + 2iu ldu 
dx 2 \dx ) 2u{u — 1) it — 1 x x dx' 
characterized by the following behavior at x = 0: 

u(x) = l + 2ix- 2n + 2 \ 2 + 0(x 3 ), a;->0, (169) 

7T 

and put 

v (x) = y/u(2x) = 1 + 2ix - 2lT + 4l x 2 + O(o: 3 ). (170) 

7T 
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Then 



= XX • ( m ) 



2iv(x) — v'(x) 
v 2 (x) — 1 

Alternatively, one can use the Hirota-Jimbo-Miwa-Okomoto a-form of Painleve V 
(see [29]), 

t *l) = -16 (a-x^-U^.Y)[a-x^), (172) 




da; 2 / \ dx A \dx J 

and choose the (regular for all positive x) solution cr(x) satisfying the initial conditions 

2 4 

a(x) = x -x 2 + 0(x 3 ), x^O. (173) 

The relation of a(x) to the function is given by the formula (see [12_Jjj]), 

It also worth noticing that the function o~(x), similar to the function £(z), can be 
determined via the solution m(z; x) of the Riemann-Hilbert problem (11651) via the 
equation (see (4.55) of [12J ) 

a(x) = —2ix(mi(x))n. (175) 

The central role of the function £(x) in the analysis of the determinants D±(x) is based 
on the following important fact. 
Define (see p] (4.38)) 

Q ± {x):=e{x)±i'{x). (176) 

Then, 

g ± (x) = -2^(log J D ± (x)). (177) 

This equation is proved in [T2] (see equation (4.125) of that work) using Dyson's results 
[T6] concerning the spectral analysis of the 1-D Schrodinger operators with the potentials 
determined by the second logarithmic derivative of the determinants D±(x). In the 
Appendix, we give an alternative derivation of ( 11771) based solely on the Riemann- 
Hilbert problem (11651) . 

Equations (11 77ft are companion equations to the equation 

e(x) = -^\o g P^, (178) 
dx A 

which in turn follows from the relation 

a(x) =xL\ og p(-). (179) 
dx 

This is one of the key formulas concerning the sine-kenel determinant P". It was first 
discovered by Jimbo, Miwa, Mori, and Sato in [30] ■ In [12] it was re-derived using the 



II In [E], the symbol 6(x) is used instead of cr(x), the symbol y(x) is used instead of v(x), and the 
symbol w(x) instead of u(x). 
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Riemann-Hilbert problem (I165p (see also Appendix). One more derivation of (j!79j) was 
obtained earlier by Tracy and Widom in |43j. 
We shall also need the important formula 

P<d = D+D-, (180) 

whose proof via the general operator technique is given in [38]. An alternative proof of 
( 11801) via Riemann-Hilbert techniques is presented in [12], page 206. 

We now proceed to establish the above mentioned evaluation of \ogD±(x) and the 
asymptotics of these functions in terms of the integrals of the Painleve V transcendents, 
i.e. in terms of the function £(z). 

Lemma 5.2. 



log .D± (x) = - 



X 

T 



log x log 2 3 , 



+ 



Proof. From [T2], page 206 we have 



x=0 



and 



and S' 



_ p0») 
dx 

d / D, 

dx" log ^ 



dx 



(logD+D- 



24 



2 

x=0 71 



ay)dy + o{\). (181) 



(182) 



x=0 



dx 



log A 



x=0 



2 

7i 



— and -^-(log£)_ 

7r dx 



0. 



x=0 



Integrating (11771) twice and using (11841) and D±(0) = 1 gives 

\ogD + {x)=-- / ?( s )dsdy-- / ^(s)dsdy--x, 
2 Jo Jo 2 Jo Jo 71 

y rx ry i rx ry 

logD_(ar) = - - / / e(s)dsdy + - / ?(s)dsdy. 
z Jo Jo 1 Jo Jo 

In view of (I178p . the integral involving £ 2 (s) can be expressed in terms of the determinant 

P^ x \ Indeed, taking into account the first equation in (11821) and the equation P = 1, 

we derive from (1178ft that 



(183) 

(184) 

(185) 
(186) 



'0 jo 
Simultaneously, from 

-x py 



\ 2 (s)dsdy+-x = -\ogP {x \ 



;i87) 



and (11711) it follows that £(0) = 2/tt, and hence 



Jo 



e(s)dsdy 



71 



% Equations (|182p follow also from the direct small- x expansions of the Fredholm determinants D±(x), 
which can be easily obtained with the help of the identity 



log£>± = trace log (l - K^ ) 
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Combining (I185p . (11861) . (11871) . and f)188p we arrive at the following representations for 
the logarithms of the determinants D±(x): 

\ogD ± {x) = l - logPW t(y)dy. (189) 

The asymptotic expansion 

logpW = ~ - I bgx + ^ + 3C'(-1) + o(l) as x -> +oo (190) 

was conjectured by Dyson [16J and proven by three different methods by Krasovsky |35j, 
Ehrhardt |Hj, and Deift, Its, Krasovsky, and Zhou [10J. Substituting ffT90|) into ffT89l) 
we obtain fl!81[> and complete the proof of the Lemma. □ 

Now we compute f£ £(y)dy. 
Lemma 5.3. We have 

f x loe 2 

U(0,x):= £{y)dy = x 1- + o(l) as x -> +oo. (191) 

Jo 2 

Proof. If m(z) satisfies the Riemann-Hilbert problem (11651) then := ^(^e 12 ^ " 3 

satisfies (see (4.11) and (4.30) in [15]) 



^ = (^3 + ^)^. (192) 

ax 

The function ip(z) solves the Riemann-Hilbert problem ((4.2) in [12J) 
' ip is analytic for z ^ J 

^ = ^ ( ^ 2 ) on 7 (193) 

?/>(;z)e- ia;ZCT3 = / + O Qj as ^ -> oo. 
Set 

i2(a;) :=<0+(O;a;). (194) 
Since R(x) satisfies the differential equation 

£ = (195) 

we have 

for some constant matrix C. Note C = -R(O), which can be computed exactly. From 
Lemma 4.5 in 1121. 
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where the principle branch of log (^4) is chosen. Since lim log 
the limit is taken from the upper half-plane, we have 



z + 1 
z- 1 



C = R(0) = lim ip(z;x = 0) 

2^0 + 



1 1 

2 V -1 3 



Hence 



R(x) 



1 / e -u(o,x) 2e c/(0 ' a:) - e~ u (°' x>> 

2 I - e -u(o,x) 2 e c/ (°' x ) + e - u (°> x ) 
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—in, where 
(198) 

(199) 



Now lim R(x) is analyzed via the nonlinear steepest-descent method for Riemann- 

x— >+oo 

Hilbert problems as in [12J. Define g(z) := \J z 2 — 1 with branch cut J and g(z) ~ z as 
z — > oo. Then 

/(z;x) := ^(z; x)e- ix9(2),J3 
satisfies the Riemann- Hilbert problem 



(200) 



' / is analytic for z ^ J 
/+ = /- 



1 

-1 2e 2ixff+ ^ 



on J 



(201) 



f( z ) =I + o[-\ as : 



oo. 



Now since 3(g+(z)) > on J, the jump matrix in (12011) decays to 



1 
-1 



as 



x — > +oo. Care must be taken because the decay is not uniform in x near ±1. 
Nevertheless, it is shown in [12] that lim /+(0; x) = /?°(0), where f°°(z) is independent 

x— >+oo 

of x and satisfies the Riemann-Hilbert problem 
' f°° is analytic for z ^ J 

1 
-1 



on J 



(202) 



= J + O ( - ] as ^ 



oo. 



Define 



:= 



1/4 



(203) 



with branch cut on J and so that (3{z) = l+0(-) as z — > oo. Then the Riemann-Hilbert 
problem (12021) is solved by 



K/3 + /3- 1 ) ^-/T 1 ) 
-K/3-/3- 1 ) K/3 + /3- 1 ) 



(204) 
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Using g+(0) = i and /3+(0) = e 17r//4 gives 

yft( 1 l\ ( 205 ) 



lim R{x)e XC73 = lim lim ^(z; x) e - i:r9{2)CT3 

;^+oo x— >+oo >0+ 



= lim / + (0) = /«(0) 

a;-»+oo Z \ — 1 1 

The (11) entry of this equation shows 

e -u(o,x)+x = ^2 + o(l) as x -> +oo, (206) 

and so 

r* 1 

/ £(y)dy = x log 2 + 2vrim + o(l) as x -> +oo (207) 

Jo 2 

for some m G Z. Since the left-hand side of (11811) is real, m = 0, which establishes 
(ELD). □ 

Taken together Lemmas 15.21 and 15.31 immediately prove Theorem 15.11 Simultane- 
ously, we have obtained the following Painleve V analog of the Painleve II total-integrals 
theorems of Sections [2] and [3j 

Theorem 5.4. [A special fifth Painleve transcendent] Suppose that u(x) is a 
solution to the Painleve V equation $168\) characterized by the Cauchy condition $1691) . 
and let £(x) be the function defined by u(x) according to equation $171}) . Then 

"OO 1 

(l-£(y))dy = -log2. (208) 

The function £(x) can be alternatively defined by equation $17 J^) in terms of the solution 
c(x) of the a -version $172$ ) of the fifth Painleve equation characterized by the initial 
condition $173). 



6. The first integrals of the mKdV equation 

The evaluation of the total integrals of the global solutions of Painleve equations 
performed in the previous sections was based on the analysis of the solution \&(A;a;) of 
the relevant Riemann-Hilbert problems at the point A = 0. One can wonder then what 
would come (if anything) from the investigation of the higher terms of the expansion 
of the \l/-function at A = oo. It turns out that if we look at these terms then, instead 
of the total integrals of the Painleve functions themselves, we will be able to evaluate 
explicitly the (properly regularized) total integrals of certain polynomials of u and its 
derivatives that play a central role in the theory of the modified Korteweg-de Vries 
(mKdV) equation 

u t — 6u 2 u x + u xxx = 0. (209) 
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We remind the reader (see [2]) that the second Painleve transcendents provide this 
equation with the important class of self-similar solutions. Indeed, if u(x) is a solution 
of the Painleve equation (TjQ) then the formula 

U 0M) = 7^T7T n ( T^YTa ) ' ( 21 °) 



(3£)V3 V( 3 *) 1/3 
gives a solution of the mKdV equation (12091) . 

A fundamental fact about equation (I209p is that it defines (for more detail see 
e.g. PQ, [IH]) on the proper functional spaces, e.g. on the Schwartz space, an infinite- 
dimensional completely integrable Hamiltonian system that possesses an infinite number 
of independent and commuting first integrals of the form 

POO 

I n = a 2 ndx, n = 1,2,3... (211) 



Here, each conserved density cik is a polynomial of u and its derivatives up to the order 
k that can be found explicitly via the following recurrence relations: 

a = - l -u\ (212) 

1 . . 

oti = ~-uu x , (213) 

a 2 = \ {uuxx ~ m 4 ) , (214) 




a k+l = 7y \ — a k - + 2^ ai0irn J ' k - ( ^' 

l,m>0; l+m=k-l / 

It can be observed that all the a's with odd subscripts are total derivatives (cf. (I213P ). 
hence the appearance of only a 2n in the description of the nontrivial first integrals (I2lip . 

Suppose now that u(x) is a global (for real x) solution of the Painleve equation ([T]). 
Then each a& can be transformed to a polynomial of u, u x , and x, 

ak = a k (u,u Xl x). 

Our aim in this section is to evaluate the properly regularized total integrals of 
ak{u,u x ,x). Obviously, we need to concentrate on the a's with even subscripts only. 
The remarkable fact is that, when calculated for the Painleve functions, the a^n become 
total derivatives (of certain polynomials of u, u x , and x) as well. This is well known in 
modern Painleve theory (see, e.g., and [37]). We shall now outline the procedure of 
finding the relevant antiderivatives. To this end let us recall the origin of the recurrence 
system fl2T2l)-fl2T5l). 

Even before their emergence in soliton theory, the densities ctk(x) were very well 
known in the scattering theory of the Dirac equation ([9]) with a rapidly decaying 
potential u(x) (see e.g. [36] and earlier references therein). The densities ak(x) appear 
in the so-called trace formulae which equate integrals ( I211I) with the moments of the 
logarithm of the absolute value of the transmission coefficient associated with the 
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potential u(x). We will discuss the trace formulae in more detail in Section What is 
important for us in this section is the main ingredient of the trace formulae derivation, 
i.e. the Riccati equation associated with the Dirac equation ([91). The Riccati equation 
appears after one transforms the first order matrix differential equation (Q to the second 
order scalar differential equation for the entry ^/n(A;x), 



*n )M! - —^ii,, + ( A - iA— - u ) * u = 0. (216) 
u V u J 

If we now write the function \t / n(A; x) in the form 

*u(A; x) = exp (~|a 3 - iXx - L(A; x)j (217) 

and put 

a(X;x) = L x (X;x), (218) 

then this substitution indeed brings equation (12161) to the following Riccati type 
differential equation for a (A; x)\ 

^ -2iAa-a 2 - — a + u 2 = 0. (219) 
ax u 

Assume now that the function a(A; x) admits the differentiable asymptotic expansion 
«( A ^)~E^T as A ^°°- ( 22 °) 

k=0 

This is certainly true in both cases of our interest, i.e. in the case of the Schwartz 
function u(x) and in the case of the Painleve function u(x). The recurrence system 
f l212l) - fl215p appears now as a result of the substitution of the series (12201) into Riccati 
equation (I219p . 

Let us now expand the function L(A; x) in the neighborhood of A = oo, 

L ( A ^)~E^' as A ^°°- (221) 

k=0 

We have that 

a k {x) = ^-L k {x). (222) 
ax 

A principal point now is that, in the case of the Painleve function u(x), all the coefficients 
L k (x) are polynomials of u, u x , and x. This important fact follows from the possibility, 
in the Painleve case, of using the first equation of the Lax pair ([8])- (Q, in addition to 
the second one, in order to study the asymptotic series (122 ip . 

Technically, it is more convenient to start with the asymptotic series for the whole 
matrix function \l/(A;x), 

$(A;x)e e ( A ^ 3 + as A - oo. (223) 

k=l 

The existence and differentiability of the series follows from the general properties of 
the Riemann-Hilbert problem (see e.g. [5], see also [22\). Moreover, the entries of 
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(225) 



the matrix coefficients m k are polynomials of u = 2 (mi) 12, u x , and x. The recurrence 
procedure that allows one to evaluate these polynomials is the following. 

As in [31] (see also Chapter 1 of [22]), we rewrite the formal series from (12231) as 

fc=i \ fc=i / \/c=i / 

where all the matrices A k and F k are assumed to be diagonal and diagonal-free, 

respectively. Then, from the differential equation (jSj) we easily get0 

u u x 1 / SN 

F 1 = -<ri, F 2 = — — er 2 , F 3 = -- [xu + u ) a x , 

F * = TZ i u + xu * + m2 ^) °"2, Ai = ^(m 4 + ^ 2 - u 2 x )a 3 , 
lb 2 

and the recurrence relation for the rest of the coefficients, 
4i[a 3 , F k+3 ] - kA k = kF k - (2u x Ox - 2iu 2 cr 3 )F k+1 

-Aua 2 F k+2 + 2J m ^mFu k > 1. 

i,m>l; l+m=k 

Note that taking the diagonal part of the last equation we determine for > 1 while 
the off-diagonal part yields F k for k > 4. The coefficients m^ of the original series are 
determined, once again by recurrence, via the identity (1224ft . Indeed, 

11 1 

mi = F x + Ai = -o x + -(u 4 + xu 2 - u 2 x )a 3 , (227) 

m 2 = F 2 + A 2 + \k\ + F1A1 

1 , 1 (228) 

= - (u 2 - (m 4 + xu 2 - it 2 ,) 2 ) 7 — - (u x . - m(u 4 + xu 2 - m 2 )) cr 2 , 
o 4 

and so on. 

Let us now come back to the series (12211) . The coefficients L k (x), which we have 
been after, are recursively determined from the already known m k (x) via the formal 
identity 



(226) 



exp 




E^H. (229) 



fc=l 

It follows then that all coefficients L k (x) are indeed polynomials of u, u x , and x. In 
particular, we have that 

L = -(mx)xi = -^u 4 + xu 2 -u 2 x ), (230) 

+ It is worth noticing that the polynomial 

H = i(A 1 )xi = ^(u 2 x -xu 2 -u 4 ) 

is the Hamiltonian for Painleve II equation ([T]) with respect to the usual choice of the canonical variables: 
q = it, p = u x - For more on the Hamiltonian aspects of the theory of Painleve equations, which we feel 
should have a strong relation to the topic of this paper, we refer the reader to the papers [39], [31], and 

M. 
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Li = \l\ - (m 2 ) n = - l -u\ (231) 

L 2 = - l -L 2 + L Li - (m 3 )n. (232) 
6 

Equation (12221) tells us that the polynomials L k (u,u x ,x) defined by (I229j) are the 
antiderivatives of the polynomials ctk(u,u x ,x) defined by (12 121) -( 12151) . In fact, only 
half of these relations - the ones corresponding to even /c's - are of interest; whereas 
the ones that correspond to odd /c's are just identities (cf. (12311) and fl 2 1 3 1) ) . Hence, 
the polynomials L 2n (u,u x ,x) are exactly the antiderivatives we have been looking for. 
We are ready now to proceed with the evaluation of the regularized total integrals of 
a 2n (u,u x ,x). 

Let u be either a purely real Ablowitz-Segur or Hastings-McLeod solution. Then 
we can integrate (12221) from x to +oo and obtain the relations 

r+oo 

I a 2n {u,u y ,y)dy = -L 2n (u,u x ,x), n = 0,1,2,... (233) 

J X 

In particular, the first relation reads3 

u 2 (y)dy = u\- xu 2 - u A . (234) 

Suppose now that u(x) is the Ablowitz-Segur solution ff211 - fl24l) . Then to regularize 
the above integral at x = —oo we need to subtract the term — /3/yjyj from u 2 (y). 
Simultaneously, the right hand side of (12341) satisfies the estimates 

u 2 — xu 2 — m 4 = — 2(3\x\ 1 ^ 2 + o(l) as x — > — oo. 

Therefore, for any c G 1R we have that 

£°° u 2 (y)dy + £ (u 2 (y) + -j= j dy - 2/3sgn(c)|c| 1 / 2 - 2p\x\ 1 ' 2 



2 2 4 

u x — XU — It 



(235) 



= -2(3\x\ l/2 + o(l) as a;^-oo. 
Hence, we obtain the following total integral formula: 

-+oo 



J^u 2 (y)dy + J^ L 2 ( 2/ ) + ^=jd 2 / = 2/?sgn(c)|c| 1 / 2 . (236) 



\y\ 

In the case of the Hastings-McLeod solution, we need to subtract off the term — y/2 in 
order to make the integral convergent at x = — oo. The resulting total integral relation 
reads 

£°°u 2 (y)dy + £ (u\y) + lj d y=* (237) 
* Of course, equation (I234|) can be checked by direct differentiation. 
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for any cel. Similar analysis can be performed with equation (I233p for any n, and it 
yields the total integral relation of the form 

J a 2n (u,u y ,y)dy + J fa 2n (u,u y ,y) - ~ = ~ F n(c) (238) 



for any c < 0. Here, the function F n (x) is uniquely defined by the asymptotic relatior 
L 2n (u(x),u x (x),x) — F n (x) + o(l) as x — ► — oo. (239) 

In particular, 

. 23(— x) 1 / 2 for the Ablowitz-Segur solution 
F {x) ={ ' (240) 

—x /4 for the Hastings-McLeod solution. 

In order to explicitly write the regularizing function F n (x) for large values of the 
number n one needs to know more terms in the asymptotics of the solution u(x). In the 
case of the Ablowitz-Segur and Hastings-McLeod solutions these terms can be relatively 
easily obtained from the substitution of a-priori ansatzs (whose existence is vouched for 
by the Riemann-Hilbert analysis) into the Painleve equation (CQ) (see [ll])). We also note 
that the total integral formulae, which are similar to (I236p -f l238|) . can be obtained for 
the case of the purely imaginary solutions u(x) as well. In the generic purely imaginary 
case one needs the regularization at x — +oo as well, and to determine the higher terms 
of the relevant asymptotic expansions is now a serious technical problem for large values 
of the number n. 

There is an interesting feature in which equations f!236[) - (I238j) differ from equations 
( !28l . (j4Tjl . ( l66l) . and ( |75l) describing the total integrals of the function u itself. Let 
us combine in all these equations the integral terms with the terms generated by the 
regularization procedure and use the symbol 

/■+0O 

v.p 



to denote this combination. That is, we put 



+oo 



v.p I a 2n (u,u y ,y)dy 

30 [ c ( dFJy)\ ^ 241 ' ) 
a 2n (u, u y , y)dy + J [a 2n (u, u y , y) -j^— J dy + F n (c) 

for the purely real Ablowitz-Segur and Hastings-McLeod solutions, 

r+oo 

v.p / u(y)dy 



c 



+ °° f c ( l\y\\ 

u(y)dy+ j \u(y) -\sx\J —\dy + is 1 —c\c\ l/2 



(242) 



(t In the derivation of (|238[) we need the differentiability of the estimate (|239j) which can be shown to 
be a consequence of the differentiability of the basic asymptotics for the solution u(x). 
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for the Hastings-McLeod solution, and 



+00 



v.p / u(y)dy := / u(y)dy 



+ / U(y) - lCrJ^r + 10" -f- - 1(7 — C^ 2 + UT-j- log C 



(243) 

2 % y a 3 4 

for the generic purely imaginary solution. The total integrals of the Ablowitz-Segur 
solutions do not need any regularization. Thus we have 

/+00 r+00 
u(y)dy := / u(y)dy (244) 
•00 J —00 

for the Ablowitz-Segur solutions. The point we want to make now is that while the 
regularized total integrals of the solutions themselves are nontrivial quantities depending 
on the solution integrated, the regularized total integrals of the densities a 2 n{u, u x , x) 
are all identically zerc tffl . The explanation of this phenomenon lies in the fact that 
all the polynomials a2 n {u,u x ,x) that have been calculated for the Painleve function 
u(x) become the total derivatives of other polynomials - the polynomials L k (u,u x , x). 
Therefore, the evaluation of the total integrals of a n (u, u x , x) becomes rather trivial due 
to (very nontrivial!) fact that we know the global asymptotics of the Painleve functions. 
In the case of the total integral of the function u itself the situation is much different. 
The antiderivative of u is actually given by equation ([TO]) and it is not a polynomial 
in u and u x . The antiderivative is explicit, but it is given in terms of the solution of 
the associated Lax pair. Therefore, in order to evaluate the integral of u the knowledge 
of its asymptotics is not enough. We have to know the asymptotics of the associated 
\I/-function. The latter we extract from the asymptotic analysis of the Riemann-Hilbert 
problem associated with the Painleve equatiorjjjfll]). 

Remark 1. The vanishing of all regularized total integrals of ac2 n (u,u x ,x), which has 
been established in this section, perhaps can be also derived using the meromorphicity of 
the second Painleve functions and the fact that the densities ct2n('u>(x), u x (x), x) all have 
zero residues at the poles of u{x) (Treves' type theorem fj^j )- 

Remark 2. As we have already seen, the evaluation of the total integral of the 
function u(x) is a more difficult task than the evaluation of the total integrals of 
a2n(u(x), u x (x), x), e.g. the evaluation of the u 2 (x). Even more difficult, though still 
possible (JEEl, WH > IIW> \M)> ^ s evaluation of the total integrals of the combination 
xu 2 [x) which, in case ofu{x) being the Hastings-McLeod solution, appears in connection 
with the analysis of the Tracy-Widom distribution functions . The integral mentioned is 
neither one of the polynomials «2n(w, u x , x) nor it can be extracted from the behavior of 

ff The particular statement that v.p J_ 00 u 2 (x)dx — was first pointed out to the fourth author by 
J. B. McLeod; it plays an important role in the analysis of the double scaling limit in the Hermitian 
matrix models with quartic potential - see [5] . 

f The evaluation of the integrals of ctk(u,u x ,x) also needs the asymptotic solution of the Riemann- 
Hilbert problem. Indeed, the global asymptotics of the function u(x) necessary for this evaluation are 
obtained from the solution of the Riemann-Hilbert problem. 
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the ^-function at X = 0. Evaluation of this integral involves an extra discretization of 
the original Painleve equation and an extra Riemann-Hilbert analysis of certain Toeplitz 
(W) or Hankel (Jffl$) determinants. The corresponding answer reads 



yu 2 {y)dy + 



[ (vAv) + T + ir) dy 

J-oo\ 2 8yJ (245) 

= — + - log c h —2- + C'(-l). 

6 8 6 11 8 24 s v ; 

A similar "most difficult" integral in the case of the solution cr(x) of the Painleve V 

equation we dealt with in Section [3] follows from the formulae \179) and U90\) . and it 

reads 

y> -dy + l-W+y+ )dy 



y J c V v % y (24g) 

= -2-4log|c| + — + 3C(-1). 

We bring the reader's attention to the appearance of the Riemann zeta-function in both 
equations. 



7. Trace formulae 



Let u(x) be a purely real Ablowitz-Segur solution of the Painleve equation (0Q), and 
let L(X;x) be as in (12171) from Section From the asymptotic analysis performed in 
Section [3] the following estimate for L(X;x) follows (cf. (1521) . (1541 ): 

z + 1/2 „/ 1 



L (z(— x) 1//2 ; x) =— //log' 



+ 



(247) 



cLS X 



z-l/2 ' ~ V(-^) 3/4 
-oo, uniformly for \z\ > 1. Here, we remind the reader that 

/i = -^rlog(l-N 2 ), -Kisi<l. (248) 



Expanding both sides of (12471) over the negative powers of z we arrive at the relations 

T^ L -^ X) = 2^2^m) log t 1 ~ |si|2 > + {j^w) • (249) 

or, recalling (I233|) . 

i r +o ° , . i 



\ 2n + i 

—X) 2 



a 2n (u,Uy,y)dy 



27ri2 2ra (2n 



' -log (l-| Sl | 2 )+0 / ' 



z) 3/4 / (250) 



as x 



-oo, n = 0, 1, 2, ... 



The left-hand side of (I250p suggests yet another way (comparing to the one used in 
Section [6]) to regularize the total integrals of a 2n - Namely, we can put 



a 2n {u,u y ,y)dy := lim / 

-oo x ^-°° {-X)— Jx 



a 2n (u,u y ,y)dy 



.(251) 
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From (12501) it follows then that 



r+oo ^ 1 

reg J ^ a 2n {u,u y ,y)dy = ~^~ 22n(2n + ^ log (l - \ Sl \ 2 ) 



(252) 



n = 0, 1,2, 



There is a striking similarity of relations (I252p with the classical trace formulae of 
the theory of the Dirac operator (Q with the potential u{x) belonging to the Schwartz 
class. Indeed, the Dirac operator trace formulae are 

30 | /"OO 

a 2n (u,u y ,u yy ,...)dy=-— log (l - |r(A)| 2 ) X 2n d\, 

> ^TTl J -oo [ZOO ) 

71 = 0,1,2,..., 

where the reflection coefficient r(A) is defined via the following relation (see e.g. |19j ; 
see also [TT]). 



$I 1 (A;x)$ + (A;x) = I 1 J^ 1 " ^ j , A G R, (254) 

where, in turn, <3>(A; x) is a unique solution of the Dirac equation (jUJ) satisfying the 
conditions 

$(A; x) is analytic in C \ R (255) 

and 

$(A; a;)e iAx,T3 -> J as A -> oo. (256) 

In other words, r(A) is defined through the jump matrix of the Riemann-Hilbert problem 
corresponding to the Dirac operator in the formalism of the inverse scattering problem 
(for more detail see again [19]). Equations ( 12521) can be formally obtained from the 
trace formulae f!253j) by putting in the latter 

log(l-|r(A)| 2 )=log(l-| Sl | 2 )^Q-A 2 ^), 

where 9{s) is the Heaviside step function, and replacing simultaneously in the left hand 
side the symbol j by the symbol reg j . We also note that the interval [—1/2, 1/2] plays 
a central role in the asymptotic analysis of the Painleve Riemann-Hilbert problem Q 
as x — > — oo - see Section |3j 

The relation between the formulae (12521) and ( 12531) can be made less formal if one 
observes that under the restrictions (l2Tj) on the monodromy data corresponding to the 
purely real Ablowitz-Segur case one can define the solution $(A; x) of the associated 
Dirac operator (Q that would almost have the properties (I255P and (J256]). Indeed, if we 
put 

^ 2 (A;x)ef A3<T3 for 3A > 0, 
$(\;x) = { (257) 
^ 6 (A;x)e- A a& for 3A < 0, 
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then we would have that 

$(A; x) is analytic in C \ R (258) 

and 

$(A; x)e iXxa3 — ► I as A -> oo, 3A ^ (259) 
(note the last inequality!). Simultaneously, in place of (12541) we get the equation 



|2 TT- 



$: 1 (A;x)$+(A;x) = ( 1 * x ), A G R. (260) 



Hence we arrive at the almost non-formal identification, 

r(A) = si. (261) 

In conclusion, we note that the "Painleve trace formulae" (1252 j) can be also used 
to evaluate the total integrals v.p. J «2 n dx that we worked out in Section El However, 
as with the method of that Section, in order to evaluate the integrals for large values of 



/(z) = (A(z), f 2 (z)Y := (e'-e-j 



the number n we need to find higher corrections to the estimate f)250p . 
8. Appendix 

In this appendix we present an alternative to the derivation of equation (11771) in [ 1 2J . 
Instead of using the spectral results of [16] we shall rely solely on the relation between 
the sine-kernel determinant and the Riemann-Hilbert problem (11651) . This relation 
is based on the representation of the sine-kernel K^ x \z,z') in the "integrable form"[j] 
(see e.g. Section 2 of [12] ) 

K(*\z,z') = fT{z) _ 9i f\ f T (z)g(z) = } (262) 

Z Z 

with the column vector-functions f(z) and g(z) defined by the equations 

T / \ T 

t i f \t (263) 

g[z) = \gx{z), g 2 {z) j := —[e 1ZX , -e lzx j 

A first key point is that the kernel R^ x \z : z') of the resolvent, = (1 — K^) -1 — 

1 = —(1 — K^) -1 !^) has the same "integrable" structure (see [26], [33], and Section 

2 in [12J) as the one indicated in (I262p . i.e., 

R^)^ G Y\ (264, 

z z 

where the components Fj(z) and Gj(z), j = 1,2, of the column vector-functions F(z) 
and G(z) are defined by the relational 

Fj = (l - KW) -1 /,, G s = (l - (KW) T ) _1 9 , j = 1,2. (265) 

\ The theory of "integrable integral operators" was pushed forward in [26] and built upon the ideas of 
[30] . It was further developed in [44] , [23] , [9] . Some of the important elements of the modern theory 
of integrable operators were already implicitly present in the earlier work 40J . 
§ It is shown in [12], Section 2, that the operator 1 — K( x ) is indeed invertible for all positive x. 
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In other words, the functions Fj(z) and Gj(z) are the solutions of the integral Fredholm 
equations 

FM-J 1 Kt4(z,zr)F j V)dzf = f s (z), J = 1,2, (266) 

and 

G,(z)- f K^(z',z)G J (z')dz' = g ] (z), j = 1,2, (267) 



J -i 

respectively. It is worth noticing that in our concrete example of the integrable kernel 
the following symmetry identities hold: 

F 1 (-z) = F 2 (z), F 2 {-z) = F x {z). (268) 

The second principal observation is that the vector functions F(z) and G(z) can 
be alternatively evaluated via the algebraic equations (see [26] and Section 2 of [12J ) 

F(z) = m±(z)f(z), G(z) = (m T ± (z)y l g(z), (269) 

where m(z) is the solution of the Riemann-Hilbert problem posed on the interval [—1,1] 
with the jump matrix V(z) defined by the formula 

V(z)=I-2mf(z)g T (z). (270) 

(We note that in virtue of the second equation in ( I262p . m+(z)f(z) = m-(z)f(z) and 

(rn^(z)j g{z) = ymF(z)j g(z).) By a direct calculation, using (I263p . we see at once 
that in our case 

/ Q 2ixz \ 

fW=(_ e _ 21 „ 2 J, (271) 

hence m(z) is exactly the solution of the Riemann-Hilbert problem (11651) . 

The last piece of the general theory that we will need is the inversion of equations 
( 12691) . i.e., the formula expressing m(z) in terms of F(z) (see again, e.g., Section 2 of 

my- 

r 1 d?' 

m(z) = I- F(z')g T (z')- , ^[-1,1]. (272) 

7-i z — z 

From this equation it follows, in particular, that the matrix coefficient mi{x) in the 
expansion (I166P admits a representation in the form 

mx{x) = j F(z)g T (z)dz. (273) 

Consider now the determinants D±(x) and let us try to evaluate their logarithmic 
derivatives with respect to x following the same line of arguments as presented on pages 
167-168 of [12J. We have 

d -logD±(ar) = -trace ( (l - K^V'-^K^ ] . (274) 



dx \\Jdx 
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A simple calculation shows that 

— = — ( cos x(z — z') ± cos x(z + z')) 

dx ttV J (275) 

= i(/ x (z) ± / 2 (^)) 5l (^) =F i(A(z) ± / 2 (z))ffi(«). 

On the other hand, taking into account the symmetries (I268p . the integral equations 
can be rewritten as 



(278) 



Fi{z)- K (x \z,z')F 1 (z')dz' - K^(z,-z')F 2 (z')dz' = f x {z) (276) 
Jo Jo 

and 

F 2 (z)- f K^{z,z')F 2 {z')dz' - f KM(z,-2f)F 1 V)dz' = f 2 (z). (277) 
Jo Jo 

By summing and subtracting (12761) and (12771) . we obtain the integral equations for the 
combinations \Fi(z) ± F 2 (z) j : 

(f 1 {z)±F 3 (z)) - J" K^\z,z')(F 1 {z!)±F 2 {z'))dz' 

= ( K f 1 (z)±f 2 (z)). 
From these equations we read that 

(l - Ki x) ) -1 (/! ± f 2 ) =Fi± F 2 . (279) 

Equations (12791) and (12751) imply that the operator ^1 — ^K±^ has kernel 

i(iq(z) ± F 2 {z)) gi {z') =F i(F x {z) ± F 2 {z))g 2 {z'), 
which yields the formula 
-r~ \°gD±(x) 

,i ,i (280) 

= -i J (F 1 (z)±F 2 (z))g 1 (z)dz±i J (F l {z)±F 2 {z))g 2 {z)dz. 

Taking into account the symmetry relations f!268j) one more time (and similar symmetries 
for gi(z) and g 2 (z)) we rewrite (12801) as 



Ai ogjD± ( x ) = _i j F 1 {z)g 1 {z)dz±i J F 1 



{z)g 2 {z)dz. (281) 



With the help of the identity (12731) . we transform (I28ip into the relation 

^\ogD ± (x) = -i((mi(a:))ii =F (mi(i)) u ). (282) 

To complete the proof of (11771) we only need to recall definition (I167P of the function 
£(x) and notice that (11741) and (11751) lead to the equation (cf. (4.33) of [12J) 

2i— (mi(a:))ii = £ 2 (a:). (283) 
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